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Abstract 



We calculate low-frequency noise (LFN) in a quantum point contact (QPC) 
^ I which is electrostatically defined in a 2D electron gas of a GaAs-AlGaAs het- 

\ erostructure. The conventional source of LFN in such systems are scattering 

potentials fluctuating in time acting upon injected electrons. One can dis- 
criminate between potentials of different origin - noise may be caused by 
the externally applied gate- and source-drain voltages, the motion of defects 
with internal degrees of freedom close to the channel, electrons hopping be- 
tween localized states in the doped region, etc. In the present study we 
propose a model of LFN based upon the assumption that there are many 
dynamic defects in the surrounding of a QPC. A general expression for the 
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time-dependent current-current correlation function is derived and applied to 
a QPC with quantized conductance. It is shown that the level of LFN is 
significantly different at and between the steps in a plot of the conductance 
vs. gate voltage. On the plateaus, the level of noise is found to be low and 
strongly model-dependent. At the steps, LFN is much larger and only weakly 
model-dependent. As long as the system is biased to be at a fixed position 
relative the conductance step, we find that the level of noise is independent of 
the number of conducting modes. From numerical calculations we conclude 
that the level of noise approximately obeys a power law as a function of fre- 
quency for frequencies larger than a threshold. At the steps for frequencies 
larger than the minimal transition rate for the dynamical impurities, we have 
S{uj) oc We are convinced that noise measurements will play a cru- 

cial role in the course of investigating the effect of the environment in QPCs. 
PACS numbers: 72.20.-i, 72.70.+m 



I. INTRODUCTION 

Electron transport through ballistic quantum point contacts (QPC), shown in Fig. |1|, 
has been extensively studied during the last decade (for a review see Ref. [|1|). One of the 
most interesting properties of these systems is the quantization [0-^ of the conductance as 
a function of the applied gate voltage (for a review see Ref. |^ and references therein). In 
addition, small oscillations in the / — Vsd curves have been found. The quantum channel 
is assumed to behave as a wave guide where the number of transverse modes is dependent 
on the gate voltage Vg as well as on the source-drain voltage Vsd- According to the present 
understanding (see e.g. Ref. |jl|), the steps in the conductance when measured vs. gate 
voltage are due to switching of the effective number of transverse modes in the channel. The 
small-scale oscillations in the / — I^^rf-curves are also believed to be caused by this kind of 
switching ^-0. 
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Low frequency noise in QPCs may originate from time- dependent internal as well as 
external parameters. Fluctuations in the potentials applied to the electrodes are examples 
of external sources of noise. Internal noise can for instance be generated by the motion of 
impurities within the contact region, or by a time-dependent rearrangement of charge in the 
doped region of the structure. 

Here we focus on noise of internal origin which is intrinsic for any realistic structure. It 
is well established that there is some disorder in the vicinity of small devices even if they 
are of high quality. In any disordered system, defects with internal degrees of freedom are 
present which may change their position due to interactions with a thermal bath. There will 
hence be a time- dependent impurity potential. These defects usually switch between two 
states, leading to a telegraph-like noise. Dynamical defects of this type were observed by 
several authors in the surrounding of even very high quality classical point contacts |8|-[T3|. 
The role of the defects with internal degrees of freedom in the surrounding of point contacts 



has been extensively discussed in the literature (for a review see Ref. 0] and references 



therein). Recently, new features have been studied [p!5|-[r?| which are due to the Kondo 
effect in scattering by two-level impurities. 

It is quite reasonable to believe that such defects, so-called elementary fluctuators (EFs), 
also exist in the surrounding of ballistic quantum point contacts. There are experimental 
results for QPCs which can be explained by the presence of defects with internal degrees 
of freedom. In particular, we have made an attempt |TB| to find theoretical support for 



experimentally observed correlation between telegraph-like noise and current in QPCs ^9 
An oscillating response to small source-drain voltages was found [^,^. 

The microscopic structure of the EFs is not yet completely clear. One of the proposed 
explanations for the presence of two-state defects is disorder-induced soft atomic vibrations. 
For low excitation energies the vibrations are strongly anharmonic and can be described 
as an atom or group of atoms moving in an effective double-well potential. Such entities 



are known as two-level tunneling systems (TLS) |^,23]. These are responsible for the low 



temperature properties of glassy materials. The generalization of the TLS model for higher 



excitation energies was worked out in Ref. |2^. Dynamical defects produce electric fields, 



slowly varying in time. The fields will scatter conduction electrons, thereby creating a sort of 
'noisy environment'. Another possibility for having fluctuations in the scattering potential 
is electrons hopping between adjacent impurity centers in the doped region. 

Our aim is to analyze low-frequency noise (LFN) in the current, caused by a 'noisy 
environment' of a QPC In this study we restrict ourselves to systems with a wide spacer 
between the doped region and the 2DEG where the quantum channel is defined, see Fig. |I|. 
If the spacer is wide enough (which is the case for high-mobility structures), the correlation 
length of the impurity potential may exceed the size of the QPC. Hence, the impurities in 
the doped region will contribute to the smooth static potential felt by the electrons in the 
quantum channel. A time-dependent rearrangement of the charge in the doped region leads 
to almost homogeneous fluctuations of the potential in space even if the voltages between 
the gate electrodes are kept constant. Effectively, these fluctuations manifest themselves as 
variations in time of the effective gate voltage. 

At the plateaus in conductance such fluctuations cannot change the current considerably. 
However, close to the steps transverse modes may be turned on or off. As a result, the level 
of LFN should increase significantly in the vicinities of the conductance steps. 

The conductance of ballistic channels is a strongly non-linear property. Hence one can 
expect to find rather unusual and informative properties of LFN of such contacts. Ex- 
perimental studies may serve as tools for extracting the principal interaction mechanisms 
between electrons in QPCs and their environment. Studies of this kind may also allow us 
to draw conclusions about the spatial distribution as well as the intrinsic dynamics of the 
impurities. 

This article is organized as follows. In Section ||, the current through an adiabatic 
contact is discussed which provides us with the starting point of our derivation of the level 



of noise. The model for low-frequency noise (LFN) in a QPC is formulated in Section IH. 
In this model, LFN is caused by EEs distributed randomly in the vicinity of a QPC. 

General analytical expressions for low-frequency noise as well as results from numerical 



calculations, will be given in Section [IT| . In Section |I^ we discuss our main results and draw 
conclusions. We have also submitted a derivation in Appendix ^ explaining why, under 
some circumstances, noise is independent of the number of propagating modes. 



II. MODEL 

The system under study consists of an adiabatically smooth channel connecting two 
equilibrium reservoirs. It is formed in a 2D electron gas by the split gate technique. The 
gates are assumed to provide a hard wall potential in the transverse direction (y), slowly 
(adiabatically) varying along the channel (x). Accordingly, the WKB approximation for the 
electron wave function is applicable. This wave function can be written as 

^{x, y,E) = Y^ a^Xni^, V, E), (1) 

where 



Xni^^y^E) = ^Jkn,\\iE,T00)/kn,\\iE,x)(f)n,x{y) (2) 



X exp 



J=Foo 



(3) 

kn,\\{E,x) = kpJe - en,±-'T{x,t). (4) 



Here E = eEp is the total energy of the electron — Ep being the Fermi energy in the leads 
at zero bias voltage, Vgd = — while kp = ylmEp/Ti j is the Fermi wave vector and 



kn,\\{E^x) the longitudinal wave vector along the channel. The transverse part of the wave 
function, (pn,x{y), depends parametrically on the longitudinal coordinate x; the correspond- 
ing 'transverse' energy eigenvalue is en,± and is measured in units of the Fermi energy Ep. 
Hence, 

= ImJ ' 

where d is the width of the channel, for simplicity assumed to be constant. 



The plus sign (+) corresponds to propagation from the left to the right reservoir while 
the minus sign (— ) corresponds to transmission in the opposite direction. 



The dimensionless function T{x,t) is the total potential in units of Ep felt by electrons 
moving in the electric field produced by the electrodes and the impurities in the doped 
region. We specify T as a sum T(x, t) = uo{x) + 6u{x,t) + vs{x) where uo{x) ~ eVg/Ep 
comes from the static part of the gate voltage. The fluctuating contribution from impurities 
and defects in the doped region is modeled by 6u{x, t), while v is the potential caused by the 
source-drain voltage {v = eVgd/^E-p). In order to match the Fermi levels in both reservoirs to 
the left and to the right, we must require that s(±oo) = ±1. To simplify we have considered 
the case of small source-drain voltages and put s{x) = 0. We also neglect the x-dependencies 
of the potentials Uq and 6u giving 

T{x,t) =uo + 6u{t). (6) 

For convenience we introduce a dimensionless parameter 77, 

[n/r]{uo)]'^ = en,±- uq. (7) 

The use of the parameter rj allows for a particularly simple interpretation in terms of the 
well known conductance quantization with gate voltage ||2^. If rj is an integer, say m, mode 
number m is just about to start conducting. Small fluctuations in the gate and source-drain 
voltages may then induce variations in the conductance between m — 1 and m in units of 
the fundamental conductance quantum 2e^ /h. When 77 = m + 0.5 we are exactly half way 
between two adjacent steps in the conductance when it is plotted versus gate voltage. From 
now on, uq is removed. The offset gate voltage Uq which before affected the bottom of the 
potential well in which the electrons move, is now instead included in the effective width d 
of the channel. 

Note that according to our model the fluctuating part 5u is x- independent. The reason 
for making such an assumption is that in most ballistic structures the scattering potential is 
rather soft. This is so because the impurities are located rather far from the region occupied 



by the conducting electrons (see, e.g. Ref. |26|). Since we are interested in the role of an 



average potential, we may neglect the weak spatial dependence of 5u. The noise properties 
in our model are thus entirely due to the dimensionless random quantity 5u{t). 
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The general expression for the current in the absence of intermode scattering [|l| is, 

I{Vsd,t) = ^ / dE6nF{E) Y.Tr,[E,6u{t)i (8) 

where Tn[E, 5u{t)\ is the probabihty of transmission between the two reservoirs for an elec- 
tron in mode n having energy E. According to our simplified model, which will be used in 
the our numerical calculations, we assume that all the propagating modes have the same 
transparency T at the Fermi level. Thus, 

T^Ef, 6u{t)] = f e[l - {n/T^f - 6u{t)l (9) 

where G(x) is the Heaviside step function. The transmission probability is hence dependent 
on the quantity 6u{t). The current is driven by the difference in chemical potentials of the 
Fermi-Dirac distribution functions in the reservoirs, 

6nF{E) =nF{E-EF- eVsd/2) - np^E ~ E^ + eVsd/2). 

(The total bias Vsd is here chosen to be symmetrically shared by the distributions to the left 
and to the right.) We have studied the low temperature regime (T = 0) for small source- 
drain voltages, where this difference simplifies to 5np{E) = eVsdS{E — Ep). The current will 
in this case be, 

HV,d,t) = ^^Y.Tn[EF,Su{t)]. (10) 

Since no scattering between the electronic states is taken into account, the only mecha- 
nism for changing the current is the possibility of having different numbers of propagating 
modes at different times. 

III. LOW-FREQUENCY NOISE 

A. General consideration 

We consider fluctuations in 6u{t) as a random process composed of contributions from 
many elementary fluctuators. Each of them, let's say the i-th one, is described by a nor- 
malized random variable ^j(t) changing between the two values ±1 with the corresponding 
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transition rates Tf and F". Within the contact, this EF will create a potential of strength 
for the electrons. Transitions between the two different states are assumed to be induced 
by interactions with a thermal bath. The rates F^ and F~ are in general different and 
determined by the nature of hopping and the interaction between the EF and the thermal 
bath. From the detailed balance principle we find, 

Fr/F+ = exp(-AABT), (11) 

where T is the temperature and A the energy difference between the states of the EF. For 
high enough temperatures the rates are almost equal, while at low temperatures there will 
be a significant difference between the two. Here we consider the case A <^ ksT <^ Ep, 
giving F~ = F^ = Fj. The dependence of Fj on A and T is determined by the hopping 
mechanism (see discussion in Ref. fl^). If the transitions are due to quantum mechanical 



tunneling the rates are independent of temperature and we have Fj oc A'^. For EF-phonon 
interaction, k = 3 while for EF-electron interaction, k = 1. If transitions are induced by 



activation, Fj oc exp{—W/kBT) where W is some activation energy [|T4|. 

The number of fluctuators has been assumed to be large enough for substituting fixed 
values of the parameters for each fluctuator with values obtained from continuous distribu- 
tions over a and F. We assume these distributions to be uncorrelated, V{a, F) = Pa{a)Pr(T). 
The nature of the statistical process will be completely defined by the distributions over a 
and F and the number of fluctuators (N). Whether or not the noise in practice is caused by 
internal or external parameters is actually irrelevant. The result we have obtained applies 
to all cases when electrons in the channel experience a time dependent potential of the kind 
Su{t) = Y^iLi O'i^ii't) , where the distributions of the random quantities Oj and Fj are to be 
discussed below. 

The gate potential 6u{t) (see Sec. ||) is composed of the electrostatic potential from 
several EEs in the doping layer. As an approximation we may nevertheless assume the 
EFs to be uniformly distributed in a spherical shell Vjam r < rmax (^ being the radial 
coordinate measured from the center of the channel) in the doped region. The probability 
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for finding an EF in a shell between r and r + dr is then P{r)dr oc r'^dr. If we subtract the 
static background charge, we obtain a dipole potential, with a fluctuating sign from each 
EF (see e.g. Ref. [Q). This implies that the interaction parameter a will be proportional 
to the inverse of the cube of the separation, a oc r~^. Normalizing under the assumption 
ctmax/omin ^ 1, wc find the distribution of P{a) to be 

Paia) = (12) 

Even if the upper limit does not enter the expression for P{a), it is important to include it in 
the following way. There must be an upper limit of possible fluctuations in the gate voltage 
since the number of fiuctuators is finite. It will be given by A^Omax- For larger gate voltages 
the one- and two-point probabilities must be zero. But the approximation we have adopted 
to find these probabilities (see below) is only asymptotically (A^ oo) exact, and disregards 
this important fact. Therefore we have neglected the finite value of amax throughout the 
derivation but included it at the end by truncating the two probability distributions over 
fluctuations in gate voltages. Note that while amin is related to the maximum distance to 
the impurities, Omax is related to the minimum distance. 

The transition probabilities usually have an exponential dependence on the barrier 
parameters. Indeed, assume the EF to be described as an entity moving in a double well 
potential [E^l-E^. The logarithm of the transition probability is then proportional either 



to the barrier height (in the case of thermal activation), or to its strength (in the case of 
quantum tunneling). Since these heights have been found to have little spread among the 
different fiuctuators, the distribution of ln(F) is smooth and almost uniform. We may hence 
model this distribution to be constant in the region between Fmin and Fmax- As the result 
we obtain 

^■■(r' = wr \r . ) P 

f-'-yJ- max/ ^ mmj 

for the normalized distribution of F. 
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B. Analytical expression for noise intensity 



Noise is usually characterized by the current-current correlation function, 

S{T) = {I{t + r)I{t)),-{m)l (14) 

or by its Fourier transform with respect to r, S{uj). The symbol (■ ■ ■)t means average over 
t which is (under stationary conditions) the same as an ensemble average over the random 
process C,(t) (ergodicity). 

From expression ([1^) for the current, we obtain the following current- current correlation 
function, 

\ / n.m 

-(T„(t)),(r„(t))j. (15) 

The transmittance is T„(t) = Tn[E-p 6u{ty\. In order to simplify the notation, we from now 
on substitute 5u with u. Since the gate voltage fluctuations must become uncorrelated for 
large time differences, expression (15) may be written as 



S{r) = ( ) E [G{n, m\T) - G{n, m|oo)] , (16) 

\ / n,m 

G{n,m\r) = {Tn{t)T^{t + r)) ^ 

= J du J dvTn{u)Tmiv)P2{u,v\T). (17) 

The two-point probability of having gate voltage fluctuations u and v at times differing 
by r is P2{u,v\t). P2 must split into one-point probabilities Pi for large values of r, 

\imr-,ooP2{u,v\T) = Pi{u)Pi{v). 

The probabilities Pi and P2 may be calculated using a generating function, 

ir;v(a;,2/) = (e-^^"W-*^^(*+")) , (18) 



where N denotes the number of fiuctuators. The average is over all random quantities 
^j(t) representing elementary fiuctuator i with strength and transition rate Fj. From the 
definition of expectation values, 
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KN{x,y\T) = J e-'''''-'y''P2{u,v\T)dudv, (19) 

Kn{x, 0\0) = J e-'''''Pi{u)du. (20) 

These relations can easily be inverted to find Pi and P2 as functions of K^, 

P2iu,v\r) = J e""+*^^K^(x,i/|r)rfxd?/, (21) 



Piiu) = ^ / e^""ir^(x,0|0)dx. (22) 



(2vr)^ 
1 

2^ 

Provided the elementary fiuctuators are statistically independent, i^Ar transforms into a 
product, 

TV 

i^^(a;,y|r) = n(e^""'^'^*^''''""^'^*^"^) • (23) 
i=i ^' 

Now consider factor number i only. This is a strictly one-fiuctuator quantity, and instead 

of fixing the values of the constants and Fj we may introduce distribution functions over 

a and F. Since then we have a function of these two parameters we have to perform a 

weighting over these distributions as well as the average over the statistical variable ^, 

{K^{x,y\T)),= {K,{x,y\T))^, (24) 

where (■ ■ ■)£ means an average over parameters a, F describing the environment, while 

Ki{x,y\T) = ^e-^^<W-'?'<(*+"))^ . (25) 

The last expression depends on the parameters a and F. In order to evaluate we apply 
the Holtsmark procedure P3 ? which is applicable in the limit of many fiuctuators. 



KN{x,y\T)=e-^^^^^y\^\ 

FN{x,y\T)=N[l-{K,{x,y\T)),]. (26) 

This approximation is valid provided 1 — Ki{x,y\T) ^ 1 and ^ 1. The generating 
function Ki{x, ?/|r) for one fiuctuator and fixed values of a and F has been calculated earlier 



fsee also Refs. [29,1 
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Ki{x,y\T) = e ^'"^l [cosh(r|r|) cos (a(x + y)) 

+ sinh(r|r|)cos(a(x-?/))]. (27) 

As is clearly seen, this expression decouples for large times into a product of functions of 
X and y as it must if gate voltage fluctuations at large time differences are uncorrelated. 
Equation ( pTj) implies that Fn in Eq. ( p6D will be given by 

FM{x,y\T) = -N mx + y) + $(x - y) 

+ *(r)[$(x + ?/)-$(x-?/)]}. (28) 

The functions $ and are here defined as 

X /"Umax 

^{x) — -Pa(a)[l — cos(ax)] da 

2 -rfa , (29) 



2 v/onjin 

amax 1 — COS (ax) 



^(r) = r"""' Pr(r)e-2ri-lrfr. (30) 

Here we have made use of the distribution over fluctuator strengths (|1^). Using (p!3D , \E'(r) 
can only be calculated analytically in the asymptotic limits, 

I 2|r|rmax I I ^ -p-1 

ln{rn,ax/rn,in) ' I ' I ^ 

■^(r) = i '';(Vkir^in) r^i « « r^i . (31) 

m(rmax/rniin) ' ' ' ™™ 

0, \r\ > r^L 

The product A^amin, which in the following will be denoted as ACn/h, has the physical 
meaning of a typical shift of the QPC potential. We will show below that the typical value 
of X in the function $(x) is C^^. Hence, we replace the lower limit of integration in Eq. 
(p9|) by 0. The behavior of $(x) depends on the relation between the upper limit Omax and 
Cat. We can also replace the upper limit by infinity provided, 

Omax > C'at, or Tmin < f. (32) 

The constants rmin and f are here the minimum and average distance, respectively, between 
the QPC and the EF. Taking the derivative of $ twice, we obtain an integral representation 
of the Dirac delta function. Using its properties while integrating repeatedly we find, 
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In the opposite limiting case, 



$(x) = ^'^""" [xl, amax2: > 1. (33) 



= X^, ttm^^X < 1. (34) 

Under the conditions (|32|) F/v in Eq. ^ is found to be 

FN{x,y,T) = -CN[{l + ^{r))\x + y\ 

+ {l-^^J{r))\x-y\]. (35) 

The Fourier transform with respect to x and y needed to get P2, given by Eq. (^), is readily 
found from ( p6D and ( |35D to be a product of two Lorentzian functions. Similarly, from ( p2D we 
calculate Pi to be a simple Lorentzian function. Unfortunately, because of the slow decay of 
the Lorentzians, our approximations made so far are not sufficient for calculating the noise. 
We must also take into account that for large values of the fluctuations {u,v > Na^a.^), the 
approximation (|35D breaks down. We correct for this by introducing truncation factors T> 
for each variable representing fluctuations in the gate voltage in the probabilities, 

Pi{u) = Vi{u)C{u\2Cn) , (36) 
P2{u,v\^) = ^V2{u)V2{v)C+C-, (37) 



C^ = C 



C- = C 



f U + V 

\ 2 


C^iv(l + *)) 


f U — V 


C7v(l-^)) 


V 2 



(38) 
(39) 



^iW = j^^ , (40) 

Ciz,w) = -^^. (41) 



The function \&(t) is given by (PPD, and has the asymptotic values 1 for t ^ and for t — »• 
00. The necessary renormalization factor imposed by the truncation procedure is A/^. For P2 
this constant has to be calculated numerically, but for Pi, Mi = 2 a.Tcta.n{Namax/TT)- Note 
that such a truncation has an obvious physical meaning - the fluctuation of, say, u cannot 
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be greater than iVamax- This shift represents the case when all fluctuator contributions add 
up coherently to maximal total potential felt by the electrons. 

The final expression for the current-current correlation function will be 

^ E/ du dv 

X [P2{u,v\^{t))-P2{u,v\0)]. (42) 

The upper limits of integration are here given by tf „ = 1 — {n/rff', and P2 by ( |3^ ) assuming 
T = 1. Finally, the noise spectrum is found as the Fourier transform of S{t). 

Let us compare our result with the conventional formalism for calculating noise in elec- 
tronic systems. To be general, assume the current to be dependent on a random quantity u 
that fluctuates in time around the average value uq = {u)t- In our case it will be the effective 
gate voltage but the nature of this quantity is not of crucial importance here. As a random 
variable, u possesses certain one- and two-point probabilities Pi{u) and P2{u,v\t), similar 
to what was presented in the preceding section. Further, consider the case for which the 
current is a smooth function of u and typical fluctuations are small. Then we are allowed 
to expand the current as a function of the variable u, 

I[uit)] = liuo) + Giuo)5uit) + 0{[uit)]^}, (43) 

where the transconductance G is defined as G{uo) = dl /du\\u=u^. Note that in our case 
this should not at all be an adequate procedure - the transmission probability being one 
for propagating modes and zero for reflected ones is certainly not a different iable function 
of the gate voltage. Nevertheless we may compare the two results to see if they deviate 
significantly. The expansion Eq. (|^) directly implies that an approximate expression for 
the current-current correlation function can be written down as 

s{T) = {i{t+T)i{t)),-{im 

^ G\uo){6u{t + T)6u{t))t. (44) 

In this approach the correlation function factorizes into two parts. One describes the rigid 
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part of the system and the other is solely determined by the time-dependent random quan- 
tities. As previously, substitute 5u with u to simplify the notation. 

In our study, u is the fluctuating part of the gate voltage and G the transconductance. 
The latter can be evaluated as the derivative of the average current with respect to the gate 
voltage. From ( ]I0| ) and (^) we easily find the average current to be. 



(n^))* = ^E[^-ctan(^)+Ar/2 



(45) 



Uk = min { 1 - ( -J , Namax } , (46) 



where the sum is over propagating modes. The transconductance will then be, 

2CJ . (47) 



h—lc ■ 



Here, fcmin = [1 + vV^ ~ A^c^max] and fcmax = IvV'^ + A^o-max]; [' ' '] dcuotcs the integer part. 
r] is related to uq according to (|^). The average {u(t + T)u{t))t can be evaluated when the 
truncation is done in the variables {u±v)/2 instead of u and v. This modification is not less 



accurate than the truncation itself. Making use of (|37D but with this modified truncation 
we find. 



+ {Na 



)2 

max/ 



(48) 



l + ^(r); \l-^(r)^ 
A{x) = 1/x arctan(x), B = 4amax/7ramm, (49) 

where ^I/(t) is given by (|30D. If Na^i^^ ^ ttCn (which is also required for using the Holtsmark 
approximation [0), the obvious restriction lim^^o('w(^ + ^)uit))t > is fulfilled. Further- 
more, since {u(t))t = we must have limT-_>oo(w(t + T)u{t))t = {u(t + T))t{u(t))t = 0, where 
we have assumed that fluctuations in the effective gate voltage become uncorrelated at large 
time differences. Obviously, (H) obeys this condition. Collecting Eqs. (|4^) , (|47|) and ( ^8|) we 
may compare this conventional approach of evaluating noise to our calculation. In Fig. (^ 
we indeed find, as expected, significant deviations between the two treatments. At the sec- 
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ond conductance step, they differ approximately by a factor 25. It sliould be emphasized 
that this comparison is strongly dependent on the choice of the parameter A^amin- 

The case opposite to (|32D for which Omax > flmin, is applicable to a structure with a very 
wide spacer compared to the thickness of the doping layer. The Lorentzians C{z\w) in Eq. 
(p6|) - ( |37D are then replaced by Gaussians, 

w\/2tx 



and the quantity Cn by A^Y^ctminOmax [see Eq. (0)]. However, we do not further pursue the 
analysis of this situation. 

C. Numerical analysis 

We had to resort to a numerical calculation of the noise using a discrete cosine transform. 



In general there is a problem with discrete transforms that is called aliasing It 
occurs when the time sequence is not bandwidth limited, and shows up as an increased 
Fourier transform for high frequencies. In short, high frequency components outside the 
range where the transform can be calculated are folded into this domain. For the case of 
investigating l//-like noise this problem is severe. The contribution to the transform that 
is transferred into the region of interest would even diverge if the transform was exactly 
proportional to 1//. When the signal has a finite bandwidth this problem does not occur. 
To invent an artificial bandwidth in our case we decrease the physical parameter Fmax to be 
of the same order as the upper limit of the Fourier transform of S{t). Evidently, this will 
be allowed since the numerically calculated transform anyway will be cut in the vicinity of 
the upper limit. 

If the function of time can only be sampled in a relatively few number of points because 
of computational difficulties there is also a problem to find the transform over a wide range 
of frequencies. Basically, it is only possible to find the transform over a number of decades of 
the order [log^^Q (number of sample points) — 1]. We solve this problem in the following way. 
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The time scale is set by the upper hmit of r, Tmax- By changing its value we may successively 
map out different parts of the Fourier transform over a wide range of frequencies. For low 
frequencies {u < lOOFmin), a large value of Tmax is used and S{t) is sampled over almost the 
entire domain. 

For high frequencies though, only the most rapidly varying part of S{t) is relevant 
because the contribution from the remaining parts is averaged out when the transform is 
calculated. Hence it is sufficient to select the most rapid variations, which in our case 
are close to r = 0. Keeping the same number of sample points while Tmax is decreased, 
the resolution in time is enhanced. With an increased resolution we are able to compute 
the high frequency components. Note also that when changing Tmax we must change F^ax 
accordingly in the way mentioned above. Unfortunately, because we don't sample S{t) for 
all values of r two problems occur. Firstly, the calculated transform will have an offset. 
We correct this by adjusting the transform to the low frequency solution where Tmax is very 
large. Secondly, there will be rapid oscillations in the transform because we truncate S{t) 
with a sharp step when Tmax is decreased. This problem is solved by filtering, averaging over 
a small interval yields the correct form of the transform. In practice we convolve the Fourier 
transform using a constant filtering function with sufficiently large bandwidth (in time). 
(Note that time and frequency are exchanged here compared to the 'normal' case when we 
filter a function of time.) We repeat the calculation and decrease Tmax in sufficiently small 
steps (we used a factor of 10) every time. Finally, we put the pieces together and match 
each new section at the lowest frequency. In this way we are able to find the transform over 
a wide frequency range. Under any circumstances, the transform for u < lOOFmin is not 
affected at all by our method of calculation. A smooth matching (see Fig. §) of the pieces 
provides some numerical evidence that the method is adequate. 

The parameter {Na^i^/ Ep) describing a typical fiuctuation in the potential from all 
impurities was chosen to be equal to 5 ■ 10"'^; the maximal fiuctuation was determined by 
c^max/omin = 100. For the Holtsmark approximation to be valid a large ratio between the 



maximal and average potential fiuctuation is required, amax/A^Omin ^ 1 (see also the end of 
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subsection [IIIBD . With these parameters we may for instance have 10 impurities implying 



o-maxZ-^flmin = 10, which sccms plausible. As has aheady been mentioned, the transmission 
probabihty Tn{u) was taken to be zero for non conducting modes and unity for conducting 
ones. In this case we will find that Y^'^=iTn{u) = [i]], where [rj] denotes the integer part of 
7], which gives the number of conducting modes at zero fluctuation in gate voltage. 

A comparison of our calculation with the simple-minded conventional approach leading 
to is done in Fig. |^. Large deviations are found for this set of parameters, at the step 
Sb{uj)/ Sa{^^) ~ 25. Indeed, significant deviations are expected due to the non-differentiable 
character of the transmittance as a function of gate voltage. The step-like shape of the 
average transmittance is responsible for the strongly enhanced level at the step and the very 
low level of noise on the plateaus. 

In Fig. 1^, the logarithm of the noise is plotted in units of {2e^Vsd/h)^ versus the loga- 
rithm of dimensionless frequency w = dj/Fmin- The graphs show the noise at the second 
conductance step [r] = 2) as well as between the second and third step {rj = 2.5). Generally, 
the current fluctuations are maximal when rj is an integer. This corresponds to a step in the 
current-gate voltage curve. Since noise is just the Fourier transform of the current-current 
correlation function, we consequently expect the noise to be large for integer and small for 
non integer values of r]. This is also clearly seen in Fig. ^ 

If A^flmax < Ep/[ri] noise is actually independent of amin o,t the step. This can be under- 
stood in the following way. The two Lorentzian functions in the two point probability ( |57D 
is in this case split symmetrically into two regions of positive and negative gate voltage fluc- 
tuations. These two domains are defined by T„ being constant. The tails of the probability 
distribution only extend into the nearest plateau on both sides of the step for this value 
of ctmax- It is then irrelevant how wide the Lorentzians are (determined by Namm/ Ep), as 
long as the total weight of probability remains the same in the different domains. Only the 
probability density profile differ for different values of Najnin/ Ep. But since is constant in 
both regions, it has no effect on S. By measuring the level of noise at and between the steps 
we may infer the order of magnitude of the parameter flmin- It will be related to the ratio of 



the two noise intensities provided A^Omax < -^f/[^]- For the presented results in Fig. ^, we 
actually have A^Omax = EY/[rj\. 

Between the steps in the current-gate voltage characteristics the noise must decrease 
with amin, since the average potential fluctuations are proportional to this parameter. This 
has also been checked. 

From Fig. |^ we deduce that in all cases, the dependence is almost linear for frequencies 
larger than some typical value uo^. It is determined by the parameters r] and Namin/ Ep. For 
T] = 2 corresponding to the second step in conductance quantization, Uc ~ 2rmin irrespective 
of the value of Na^^i^/Ep when this quotient is sufficiently small (< 1 • 10^'^). If the system is 
biased between the second and third step r] = 2.5, Uc is of the order lOFmin. This threshold 
frequency has been found to slightly increase with decreasing A^amin between the steps. 
Because of computational difficulties, the form of the curves for, let's say u > lO^Fmin, 
should only be considered to be accurate within a relative error of the order 5%. The slopes 
for uj>ujc were found to be close to -0.85 (-0.93) for r] = 2.0 (r/ = 2.5), S{uj) oc i/cjO-85(o.93)_ 
This means that the noise is almost 1 / /-like. 

Experimentally, one is also interested in the integrated noise intensity per decade. For 
the curves in Fig. |^, this quantity is shown in Fig. ^ The integration is for decade number k 
performed over the interval [k — 0.5, k + 0.5]. The center frequency is here Uc, which in turn 
defines k, k = uJc/T^[^. If the noise was identically 1//, the integrated noise intensity per 
decade would be constant for all decades. In reality, the noise must be cut for frequencies 
larger than the physical value of Fmax- In realistic systems Fmax can be as large as lO^^Fmin. 
This is beyond the limits of our calculation, so we are unable to detect this physical cut-off. 

IV. DISCUSSION AND CONCLUSIONS 

We would like to stress that the calculated level of noise significantly differed from a 
conventional estimate, see Fig. (It is lower by a factor 25 for one set of parameters). The 
approximation was based upon an expansion of the current- current correlation function in 
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fluctuations of tlie gate voltage. Due to the sharp structure of the average transmittance 
{Tn)t the conventional approach leads to invalid results. The deviation from the true result 
depends on the strength of the fluctuators - the weaker they are the larger is the difference. 

According to our model the current-current correlation function S{t) (^) has the follow- 
ing properties. Under certain restrictions (see Appendix ^ it is independent of the number 
of propagating modes, [rj]. Furthermore, if the spacer is thick enough and the inequality 
Na^nax < Ep/[ri] holds, S{t) is also independent of amin at the steps in conductance quanti- 
zation. Hence, measuring the noise level at a step would be enough to find the value of Omax 
which is the interaction strength between the nearest elementary fluctuator and the QPC 

Our model (which assumes a constant transmittance T„ on every plateau) suggests that 
the level of noise is very low close to the middle of a plateau. Generally, it is more dependent 
on the maximum (rmax oc 1/amin) than the minimum (rmin oc l/flmax) distance to the EFs. 
This fact provides us with a way of determining amin almost independently of Omax- Since 
ajnin then effectively is the only parameter, the level of noise at a plateau where many modes 
are conducting should provide sufficient information to determine its value. Having obtained 
both parameters of the model one can examine its validity by comparing the predicted and 
measured dependence of the noise on the gate voltage. 

Numerical calculations indicate that there is a typical frequency, Uc, specific for each set 
of parameters and above which S{uj) oc l/cj". The exponent a was found to be slightly 
less than 1 and the threshold frequency oOc close to the minimal transition rate rmin for the 
fluctuators. Also, on the plateaus a is slightly enhanced compared to its value at the steps. 
Further, Uc increased when the minimum fluctuator strength Omin decreased. 



Experiments on l//-noise in QPC |31] are in qualitative agreement with our theoretical 
results. Indeed, the noise was found to be S^u) oc I^/u°', where /3 = 2 ± 0.4 and a = 
0.9 ± 0.1. As we also conclude, the value of a was larger on the plateaus than at the steps 
in conductance quantization with gate voltage. Of course, the level of the noise was larger 
at the steps than between them. For the highest two values of the number of conducting 
modes, 4 and 5, the noise was almost the same in agreement with our claims (see appendix 
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It was also pointed out that after some cool-downs of the sample the peaks of the 
noise as a function of the parameter t] were smeared out. As a possible cause for this effect 
we suggest that the elementary fluctuators may change their positions between subsequent 
measurements because of heating. This would affect amin as well as amax and give rise to the 
observed effect. 

Finally, we recall that our derivation of noise is done in the limit of many fluctuators, 
^ 1. At the same time, we assume that A^ <^ (?"max/''"min)^- The last inequality means 
that a significant number of EFs are situated in a close vicinity of the QPC The probability 
to have a large time-dependent fluctuation of the electric potential is large and governed by 
a Lorentzian distribution function (|37D . If the last inequality does not hold only the nearest 
EF is of importance and the distribution is a Gaussian function. If A^ is of order 1, the noise 
has a random telegraph character. 
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APPENDIX A: NOISE FOR DIFFERENT NUMBERS OF PROPAGATING 

MODES 

The magnitude of the noise will here be shown to be independent, under certain re- 
strictions, of the number of propagating modes through the contact. Let zero fluctuation 
in gate voltage correspond to a fixed position on the quantized plateaus in the current- 
gate voltage characteristics. Assume that the Lorentzians in the two-point probability 
extend over p steps in the I — Vg curve. This provides a restriction on the value of t], 
r] > p. We have then Tn{u) = 1 for all modes n < [r]] — p = q and gate voltages u 
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and V for which P2{u,v\^{t)) = P2{u,v\t) is finite. The dimensionless correlation function 
s{t) = {h/2e^VsdfS{T) (see H) can then be reduced, 

s(r) = j j dvTm{u)Tn{y)V{u,v\r) (Al) 

n,m=q+l 

where V{u,v\t) = P2{u,v\t) — P2{u,v\0). Indeed, for terms with both Ts constant, T„ = 
Tm = 1 the integral J du J dvV{u,v\T) vanishes because of the normalization condition. On 
the other hand, if only one T is constant, let's say Tm = 1 the integrals / duV{u,v\T) vanish 
because 

J duP2{u,v\T) = Pi{v) 

is r-independent. Note that T„ is dependent of the offset gate voltage which is related to 
rj. Because of the simple form of T„, being zero for reflected modes and one for propagating 
modes, we are able to make a translation in the mode index, 

Tn{u,r]) =Tn-q{u,T] - q). (A2) 

Finally, changing indices of summation we obtain. 



J du J dv 



S[T) = 

n',m'=l ' 

xT^, (u, {r]} + p)Tn' {v, {r]} + p)V{u,v\t) (A3) 

where {rj} = rj — [rj] is the fractional part of rj. Evidently, this is independent of the number 
of propagating modes, [r]]. 
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FIGURES 



FIG. 1. A quantum point contact is defined at the interface between the undoped AlGaAs 
spacer layer and GaAs by the spht gate confinement technique. Noise is caused by the motion of 
impurities in the n-doped AlGaAs layer. 

FIG. 2. Logarithm of noise, logio['S'(a; = lOOFmin)] as a function of rj. The unit of S{uj) is 
{2e^VsdT / h)^ . The parameter rj gives the operating point, i.e. where on the conductance-gate 
voltage staircase the system is centered. The integer part of r) is the number of conducting modes. 
We calculated Sa as the numerical Fourier transform of the valid expression for the noise, Eq. 
(42). Using the traditional approach, Sb was computed from Eq. (44, 47-49). The minimum 
fluctuator potential Omin was defined by A^amin/E^p = 5 x 10~^ {Ep being the Fermi energy); amax 

by amax/amin = 100. 

FIG. 3. Logarithm of noise, logjo['S'(^^)] (in units of {2e^VsdT / h)^) as a function of the loga- 
rithm of dimensionless frequency, logiQ(ti;) (uj = wT„ii^). The minimum fluctuator potential Omm 
was defined by Na-aiin/E-p = 5 x 10~^ {E-p being the Fermi energy); Omax by amax/omin = 100. The 
calculation was done for two different values of ry, 77 = 2 where mode number 2 is just about to 
become propagating, and 77 = 2.5. The numbers indicate the approximate slope k {S{u) cx ui^) at 
frequency lj = lOOFmin- 

FIG. 4. Integrated noise intensity per decade as a function of dimensionless center frequency 
Wc = Wc/Fmin • The different plots correspond to the same cases as in Fig. 3, rj = 2, and rj = 2.5. 
The integration was performed over the interval w G [wc — 0.5, Wc + 0.5]. If S{uj) cx l/w", where 
a = 1, the bars would be of equal height. 
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